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Abstract 



r^ A /3-skeleton, /3 > 1, is a planar proximity undirected graph of an Euclidean points 

r) set, where nodes are connected by an edge if their lune-based neighbourhood con- 

• tains no other points of the given set. Parameter (3 determines the size and shape 

O of the lune-based neighbourhood. A /3-skeleton of a random planar set is usually a 

disconnected graph for /3 > 2. With the increase of /3, the number of edges in the 

^-H ;5-skeleton of a random graph decreases. We show how to grow stable /3-skeletons, 

>J^ which are connected for any given value of /3 and characterise morphological trans- 

QQ formations of the skeletons governed by (3 and a degree of approximation. We spec- 

Q^ ulate how the results obtained can be applied in biology and chemistry. 

^^ Keywords: proximity graphs, P-skeletons, pattern formation, morphogenesis 

cn 



C^ 1 Introduction 



A planar graph consists of nodes which are points of the Euclidean plane and 
edges which are straight segments connecting the points. A planar proxim- 
ity graph is the planar graph where two points are connected by an edge if 
they are close in some sense. Usually a pair of points is assigned a certain 
neighbourhood, and points of the pair are connected by an edge if their neigh- 
bourhood is empty (does not contain any points of the given set). Delaunay 
triangulation [lOj, relative neighbourhood graph [16], Gabriel graph [25j, and 
spanning tree, are the classical examples of the proximity graphs. 

The ^-skeletons [20] make a unique family of the proximity graphs monotonously 
Preprint submitted to Elsevier 9 April 2013 



parameterised by /3. Two neighbouring points of a planar set are connected 
by an edge in a ;5-skeleton if a lune-shaped domain between the points con- 
tains no other points of the planar set. The size and shape of the lune is 
governed by /3. The /5-skeletons are worth to study because they are amongst 
the key representatives of the family of proximity graphs. Proximity graphs 
are applied in many fields of science and engineerings: from image process- 
ing (eps. reconstructing the shape of a two-dimensional object, given a set 
of sample points on the boundary of the object), visualisation and physi- 
cal modelling to analysis and design of communication and transport net- 
works [T]l3l[Wl[8l[MTM71[2T^^ 



A ;5-skeleton is the Gabriel graph [25] for ;5 = 1 and it is the relative neigh- 
bourhood graph for /3 = 2 [16,20]. A /3-skeleton, in general case, becomes 
disconnected for /3 > 2 and continues losing its edges with further increase 
of /5. In our previous paper [4J we demonstrated that ;5-skeletons of random 
planar sets lose edges by a power low with the rate of edge disappearance pro- 
portional to a number of points in the sets. Some /3-skeletons conserve their 
edges for any /3 as large as it could be. These are usually skeletons built on a 
regularly arranged sets of planar points, however even minuscule impurity in 
the regular arrangement of points leads to propagation of an edge loss wave 
across the otherwise stable skeleton. Can we produce connected ;5-skeletons 
for arbitrarily large values of /5? How do these skeletons look like? What are 
properties of these skeletons? How topological features of the connected /3- 
skeletons are changed with the increase of /3? We answer these questions in 
the paper. 



2 /3-skeletons 



Given a set V of planar points, for any two points p and q we define a /3- 
neighbourhood t//3(p, q) as the intersection of two discs with radius /3\p — q\/2 
centered at points ((1 — 1)]9, |g) and (|]9, (1 — |)g), /3 > 1 |[20l[T6] . see examples 
of the lunes in Fig. [T] Points p and q are connected by an edge in ;5-skeleton 
if the pair's /3-neighbourhood contains no other points from V. 

A ;5-skeleton is a graph ^^(V) = (V,E,;5), where nodes V C R^, edges E, 
and for p, g G V edge (pq) G E if t//3(j9, q) D V/{p, q} = 0. Parameterisation /5 
is monotonous: if f3i > ^2 then B^,(V) C ^^^{^f) p[6l[20] . 

A /5-skeleton is a non-planar graph for ;5 < 1, see example in Fig.[2|3. Therefore 
we consider only skeletons with (5 > 1. 




o 

(a) /3 = 1 (b) /3 = 2 (c) /3 = 10 (d) /3 = 100 

Fig. 1. Examples of lunes, /3-neighbourhoods, of two planar points, black disks, for 
various values of /3. The /3-neiglibour hoods are shaded in grey. 

3 On stability of /5-skeletons 



^-skeletons of random planar sets lose their edges by a power law when /3 
increases linearly, see details in [4j. See examples of /3-skeletons, /3 =1, 2, 3, 4, 
7, 20, 100 and edge loss curve in Fig. [2} Most /3-skeleton lose their edges with 
increase of /3 however some /5-skeleton do not. A stable /5-skeleton retains its 
edges for any value of ;5 > 1. A most obvious example of a stable ;5-skeleton 
is a skeleton built on a set of planar points arranged in a rectangular array. 

Proposition 1 Rectangular lattice is a stable /3 -skeleton. 

This is because when /3 tends to infinity the ^^-neighbourhood tends to a 
rectangular shape and becomes an intersection of two half-planes. In details, 
let Hab be an open half-plane bounded by an infinite straight line la passing 
through a, perpendicular to segment (a, b) and containing 6; and H^^ be an 
open half-plane bounded by an infinite straight line /^ perpendicular to seg- 
ment (a, 6), passing through b and containing a. Let M^^ = H^^ H H^^. When 
P becomes extremely large, tends to infinity, a ^^-neighbourhood of any two 
neighbouring points a and b tends to M^^. A /3-skeleton of planar set V is sta- 
ble if for any a, 6 G V M^^ does not contain any points from V apart of a and 
b. The rectangular /3-skeleton conserves its edges for any value of f3 (Fig. |3]). 
The rectangular lattice is stable because for any two neighbouring nodes a and 
b intersection M^^ of their half-planes fits between rows or columns of nodes 
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Fig. 2. Example of /3-skeletons of a random planar set. (a) A random planar set V 
of 500 discs, radius 2.5 each, randomly distributed in a disc radius 250. (b-h) Exam- 
ples of /3-skeletons on planar set V. (i) Example power curves of edge loss, e(n, /3) 
is a number of edges in /3-skeleton constructed on a random set with n nodes, 
1^/3^ 100, n — 350 (solid line), n — 550 (dashed line) and n — 750 (fine dashed 
line), values of j3 are incremented by 0.1. 

without covering any other nodes. 



Given /3, is it possible to generate a planar set which /5-skeleton is a connected 
graph? A method of growing such sets, and their /5-skeletons, is presented 
further. 



Fig. 3. Rectangular lattice is a stable /3-skeleton. /^-neighbourhoods, /3 ^ oc, of two 
pairs of nodes (marked by grey discs) are shown by hatched areas. 

4 Growing /5-skeletons 



A graph is connected if there is a path along edges of the graph between 
any two nodes of the graph. An indirected graph is connected if there are no 
isolated nodes. To grow a connected /5-skeleton for a given value of /3 we start 
with a single planar point po = (xo,?/o) and then introduce additional points 
one by one. When a new candidate point is introduced to a planar set we 
check if 

• the candidate point does not fell into /^-neighbourhoods of existing nodes, 
and 

• the skeleton of the planar set with the candidate point retains its connec- 
tivity, i.e. there are no isolated nodes. 

Points can be added to the planar set either in a random fashion or in a regular 
manner. We adopt a regular addition of nodes by the following procedure. 

Node Addition Procedure 



(1) r = 5,^ = 0, 5= 10 

(2) p = (xq + r cos 9^yo + r sin 9) 



(3) if B/3(V U {p}) is connected graph and (Vg G V : |p — Q'l > S) 
then V ^ V U {p} 

(4) e^e + Ae 

(5) if ^ > 360 then ^ ^ and r ^ r + Ar 

(6) go to step 2 

We use polar coordinates p = {xq + rcosO^yo + rsinO) and assume that no 
two points can he closer than 6 to each other, in all experiments 6 = 2.5. 
Position of initial point is fixed, po = (xo,2/o)- First candidate point is placed 
at distance r from po with angle ^ = 0. The angle 9 is incremented by A9. 
When angle 9 reaches 360 degrees radius r is incremented by Ar and 9 is 
assigned value 0. The iterations may continue indefinitely but in experiments 
illustrated here we stop growing skeletons when r reached 90. Further in the 
paper we sometimes address /5-skeletons grown by above procedure as simply 
/5-skeleton. 

Proposition 2 Let A9 -^ then grown 13 -skeleton is transformed from a 
hexagonal lattice^ for /3 = 1^ to an orthogonal lattice^ for /3 ^ oo. 



For /5 = 1 lune t/^(p, q) is disc with diameter pq (Fig. Mk). Thus, in principle, 
we can arrange any number of points around an initial point po. However, 
due to imposed minimal distance 5 = 2.5 between any two points, the points 
can be considered as discs. A hexagonal packing is a densest arrangement of 
identical discs. Rectangular lattice is a stable /3-skeleton, it remains connected 
for arbitrarily large /3. 



5 Dynamics of skeletons controlled by /5 



Skeletons grown on computers are never ideal, and never become rectangular 
lattices, due to impurities in their topologies introduced by increments of 9. 
Examples of /5-skeletons grown from a single point with angular increment 
A^ = 0.5 are shown in Fig. [4] A hexagonal arrangement of nodes in /5-skeleton 
for /3 = 1 is well seen (Fig. [4^). The hexagonal arrangement is gradually 
destroyed when /3 increases from 1 to 2 (Fig. |4^bc) with majority of nodes 
having three or four neighbours (Fig. [5^ 



Further increase of /5 leads to dissociation of cycles and formation of tree- 
like skeletons with domains of lattice-like arrangements (Fig. |4|-t). Sizeable 
domains of rectangular lattices are still observed at /3 = 10, e.g. domains 
located in southern, western and north-westerns parts of the graph in Fig. |4[. 
Branching of the tree is reduced with increase of /3 till skeleton is transformed 
to a cross-like structures with a single binary branching at each of four main 
branches (Fig. [ijb). This ^-induced transformation is refiected in decrease in 
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Fig. 4. Examples of /3-skeletons grown from a single point, 
r = 5, Ar = 0.5, A6i = 0.5, S = 2.5 
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Fig. 5. Average degree (a), number of nodes (b), total lengths of 
edges (c), and diameter in nodes (d) of /3-skeletons grown with parameters 
r = 5, Ar = 0.5, A9 = 0.5, 6 = 2.5. o 



average degree of the graphs' nodes, which almost stabihses around value 2 
when f3 exceeds 50 (Fig. [5^ 



These structural transformations are reflected in decrease of a total number 
of nodes (Fig. ^) and total length of edges (Fig. [5J3). Increase of /3 does not 
affect diameters of the /3-skeletons, which vary between 23 and 29 nodes for 
1 < /3 < 100 (Fig. ^ 



6 A9 influences morphologies 



Morphologies of /5-skeletons are affected not only by values of /3 but also A^. 
We can expect that with increase of A^ the skeletons converge from almost 
regular lattices to trees for /3 < 4. For example, compare skeletons grown 
for A9 = 5 (Fig. [6^-1) and A9 = 10 (Fig. [6}2i-t). Increase of A^ causes 
predominant decrease of lateral (aligned along concentric cycles centred in po) 
edges, see e.g. Fig. [6^bc and Fig. [6|iino, and decrease in branching of trees, 
see e.g. Fig. [6|kl and Fig. [6):st. New nodes are added to /3-skeletons in a cycle 
of iterations, change of 9 from to 360, thus branches of the graphs are leaned 
towards circular arrangements (Fig. |6|^lrs). 

For /5 up to 3, /5-skeletons consist of two morphologically distinctive com- 
ponents: internal core of a quasi-regular network and spider-web like halo of 
radial rays connected by lateral links. The pronounced examples of this mor- 
phological division are shown in Fig. [Gjuno and Fig. [7^b. 

Increase of A9 leads to a disappearance of lateral links and shrinking of the 
quasi-regular network core, see transition from Fig. [Hjii to Fig. [7^ to Fig. [7[. 
Skeletons grown with large angle increments A^ are characterised by a preva- 
lence of radial edges, or rays. Thus, for A^ = 30 we observe transition from 
a small spider web when /5 = 1 (Fig. [?[) to twelve rays, /5 = 2 (Fig. [^), and 
four rays, f3 > 2Q (Fig. [^), structures. Number of rays is changed from nine, 
^ = 2, to seven, ;5 = 6, to three, /5 = 9 (Fig. [7|iop). Skeletons grown with 
large angular increments, A^ = 80 and 90, Fig. [8[-q consist of three, four or 
five rays, for any value of /3. 

In situations when 360 mod A9 > lateral edges are formed between rays. 
Examples are ;5-skeletons grown with A9 = 50 (Fig. |8^b) and A9 = 50 

(Fig.iiij.) 



(a) A6' = 5, /? = 1 (b) A6' = 5, /? = 2 (c) A6i = 5, /3 = 3 (d) A6' = 5, /3 = 10 




(e) A^ = 5, ^ = 15 (f) Ae = 5, /3 = 20 (g) A9 = 5, 13 = 25 (h) A9 = 5, ^ = 30 




(i) Ae = 5, /3 = 35 (j) A0 = 5, /3 = 40 (k) A^ = 5, /3 = 45 (i) A9 = 5, P = 50 




(m) A6I = 10, /? = 1 (n) AO = 10, l3 = 2 (o) A6' = 10, /? = 3 (p) A6i = 10,/3 = 10 




(q) A^ = 10,/3 = 20 (r) A^ = 10, /? = 30 (s) A^ = 10, /3 = 40 (t) A^ = 10, /3 = 50 
Fig. 6. /3-skeletons grown with A6 = 5 (a-1) and A6 = 10 (m-t). 
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(a) A6I = 15, /3 = 1 (b) AO ^ 15, /3 ^ 2 (c) A6i = 15, /3 = 3 (d) A^ = 15, /? = 10 





(e) A^ = 15, ^ = 20 (f) A^ = 15, /3 = 30 (g) A9 = 15, /3 = 40 




(h) A9 = 15, 13 = 50 (i) A9 = 30, 13 = 1 (j) A^ = 30, ^ = 2 






(k) A9 = 30, 13 = 10 (i) AO = 30, /3 > 26 (m) A^ = 40, ^ = 1 (n) A^ = 40, ^ = 2 





(o) A^ = 40, ^ = 6 (p) A^ = 40, ^ = 9 
Fig. 7. /5-skeletons grown with A9 = 15 (a-h), AO = 30 (i-1) and AO = 40 (m-p). 
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(a) A6I = 50, /3 = 1 (b) AO ^ 50, /3 ^ 2 (c) A9 = 50, /3 = 10 (d) A^ = 50, /3 > 30 




(e) A^ = 60, /3 = 1 (f) A^ = 60, /? = 2 (g) AO^QO, /3>3 





(h) A9 = 70, /3 = 1 (i) A^ = 70, ^ = 2 (j) A^ = 70, ^ = 3 (k) A^ = 70, ^ > 




{i) A^ = 80, ^ = 1 (m) A^ = 80, ^ = 2 (n) A^ = 80, ^ = 6 (o) A9 = 80, /3 > 9 

? ? 





(p) A^ = 90, /3 = 1 (q) A^ = 90, /3 > 1 

Fig. 8. ;5-skeletons grown with AO = 50 (a-d), AO = 60 (e-g), A(9 = 70 (h-k), 
Ae = 80 (l-o) and AS = 90 (pq). 
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(c) 

Fig. 9. Slime mould P. polycephalum growing on a corn meal 2% agar gel (a), 
non- nutrient 2% agar gel (b), and a slightly wet filter paper (c). See details in ^. 

7 Possible applications in sciences 



Some aspects of a morphological dynamic of /5-skeletons, controlled by /3 
and A0, resemble substrate-induced morphological transformations in bac- 
terial [14j and myxomycetes [2j colonies. Typically, a high concentration of 
nutrients in a growth substrate leads to formation of dense quasi-uniform 
omni-directionally propagating patterns. A low concentration of nutrients in 
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a substrate leads to formation of branching tree-like structures. Two examples 
are shown in Fig. [9| 

When Plasmodium of P. polycephalum is inoculated on an agar plate with high 
concentration of nutrients (2% corn meal agar) the Plasmodium's growth-front 
propagates similarly to a circular wave. A dense network of protoplasmic tubes 
is formed inside the Plasmodium's body (Fig. [9^). Such growing pattern might 
be matched well by /3-skeletons grown by the Node Addition Procedure, given 
in Sect, g with A^ = 0.5, 1 < /5 < 3.5 (Fig. ^f), and A^ = 5, 1 < /5 < 3 
(Fig. |6^bc) and A^ = 10, 1 < ;5 < 3 (Fig. |6[nno). 



On a non-nutrient 2% agar Plasmodium forms a tree like structure (Fig. ^ 
The Plasmodium trees are alike /5-skeletons generated with A9 = 0.5, 30 < 
/3 < 50 (Fig. §kl), A^ = 5, 30 < /3 < 50 (Fig. jefi-l), A^ = 10, 20 < ^5 < 40 
(Fig. |6]:irs). A degree of branching, or 'bushiness', of protoplasmic trees de- 
creases with increase of harshness of a growth substrate. For example, Plas- 
modium cultivated on a filter paper, instead of agar gel, produces protoplasmic 
trees with very low degree of branching (Fig. [9fc). Protoplasmic trees grown in 
a harsh conditions resemble ;5-skeletons grown with A9 = 0.5, 600 < /3 < 800 
(Fig. IJst), A^ = 15, 10 < ;5 < 50 (Fig. |6|l-h). 

The procedure of growing /5-skeletons may distantly pass for, at least at a 
very abstract level, a synthesis of cyclic (Fig. [To^) and dendrimer molecules 
(Fig. [To)3c). Skeletons grown with low values of A9 and /3 may be considered 
analogous to cyclic molecules and skeletons produced with high values of A^ 
and /3 are alike dendrimer molecules. 

Whilst no direct matching between /3-skeletons and cyclic or dendrimer molecules 
can be demonstrated we calculated Randic index [28] of the /3-skeletons grown 



for /5 up to 100 (Fig. 11). The Randic index [28] is calculated as i? = ^^j Cij * 



(^=y^), where Cij is an adjacency matrix of a graph. 



The Randic index R (originally called by Milan Randic as molecular branch- 
ing index) [28j characterises relationships between structure, property and 
activity of molecular components [H]. There are proven linear relations be- 
tween the Randic index and molecular polarisability, enthalpy of formation, 
molar refraction, van der Waals areas and volumes, chromatographic reten- 
tion index [19j, cavity surface areas calculated for water solubility of alcohols 
and hydrocarbons, biological potencies of anaesthetics [18j, water solubility 
and boiling point [15] and even bio-concentration factor of hazardous chem- 
icals [30]. Estrada [12] suggested the following structural interpretation: the 
Randic index is proportional to an area of molecular accessibility, i.e. area 'ex- 
posed' to outside environment. The Randic index decreases with increase of 
j3. Exposure of /5-skeletons is proportional to /5. This is how properties of the 
molecules imitated by /3-skeletons will change when a molecule is transformed 
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Fig. 10. (a) Benzocoronene, (b) PAMAM-OS-trimethoxysilyl dendrimer. 
(c) Tetrakis(4-poly(nitriromethylene-bis(p-phenylene))phenyl)methane dendrimer 
generation 4. Chemical structures are drawn with www.chemspider.com. 



from, e.g. aromatic to dendritic. 
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Fig. 12. Evolution of grown ^-skeletons. 



8 Summary 



Given a random planar set, a /3-skeleton of the set is, in general, a disconnected 
graph for /5 > 2. We presented a procedure for growing the ;5-skeletons which 
remain connected for any, yet specified during the growth, value of /3 however 
large it is. In computational experiments we demonstrated that with increase 
of /3 and/or decrease of approximation accuracy, A0, the skeletons undergo 
a transformation from almost regular lattices or networks to branching trees 
to cross-like graphs (Fig. 12). We speculate that such evolution of the /3- 



skeletons somewhat imitates morphological transformations of myxomycetes 



16 



and bacterial colonies governed by concentration of nutrients in their growth 
substrates and transformation of molecules from aromatic to dendritic forms. 



References 



[1] Adamatzky A. Developing proximity graphs by Physarum Polycephalum: Does 
the Plasmodium follow Toussaint hierarchy? Parallel Processing Letters 19 
(2008) 105-127. 

[2] Adamatzky A. Physarum Machines (World Scientific, 2010). 

[3] Adamatzky A. (Ed.) Bioevaluation of World Transport Networks (World 
Scientific, 2012). 

[4] Adamatzky A. How /3-skeletons lose their edges. Submitted (2012). 

[5] Amenta N., Bern M., Eppstein D. The Crust and the /3-Skeleton: Combinatorial 
Curve Reconstruction. Graphical Models and Image Processing. 60 (1998) 125- 
135. 

[6] Beavon D. J. K., Brantingham P. L. and Brantingham P. J. The influence of 



street networks on the patterning of property offenses, www . popcenter . org/ 
llibrary/CrimePrevent ioii/Volume_02/06beavoii . pdf , 

[7] Billiot J. M., Corset F., Fontenas E. Continuum percolation in the relative 



neighbourhood graph. larXiv: 1004.5292 



[8] Dale M. R. T. Spatial Analysis in Plant Ecology (Cambridge University Press, 
2000). 

[9] Dale M. R. T., Dixon P., Fortin M.-J., Legendre P., Myers D. E. and 
Rosenberg M. S. Conceptual and mathematical relationships among methods 
for spatial analysis. Ecography 25 (2002) 558-577. 

[10] Delaunay B. Sur la sphere vide, Izvestia Akademii Nauk SSSR, Otdelenie 
Matematicheskikh i Estestvennykh Nauk, 7 (1934) 793-800. 

[11] Estrada E. Generalization of topological indices. Chem. Phys. Lett. 336 (2001) 

248-252. 

[12] Estrada E. The structural interpretation of the Randic index. Internet 
Electronic Journal of Molecular Design 1 (2002) 360-366. 

[13] Gabriel K. R. and Sokal R. R. A new statistical approach to geographic variation 
analysis. Systematic Zoology 18 (1969) 259-270. 

[14] Golding I., Kozlovsky Y., Cohen I., Ben- Jacob E. Studies of bacterial branching 
growth using reactiondiffusion models for colonial development Physica A 260 
(1998) 510-554. 



17 



[15] Hall L.H., Kier L.B. , Murray W.J. Molecular connectivity. II. Relationship to 
water solubility and boiling point. J. Pharm. Sci. 64 (1975) 1974-1977. 

[16] Jaromczyk J. W. and G. T. Toussaint, Relative neighborhood graphs and their 
relatives. Proc. IEEE 80 (1992) 1502-1517. 

[17] Jombart T., Devillard S., Dufour A.-B., Pontier D. Revelaing cryptic spatial 
patterns in genetic variability by a new multivariate method. Heredity 101 
(2008) 92-103. 

[18] Kier L.B., Hall L.H., Murray W.J., Randic M., Molecular connectivity. I. 
Relationship to nonspecific local anesthesia, J. Pharm. Sci. 64 (1975) 1971- 
1974. 

[19] Kier, L. B.; Hall, L. H. Molecular Connectivity in Chemistry and Drug Research. 
Academic Press, 1976. 

[20] Kirkpatrick D.G. and Radke J.D. A framework for computational morphology. 
In: Toussaint G. T., Ed., Computational Geometry (North-Holland, 1985) 217- 
248. 

[21] Legendre P. and Fortin M.-J. Spatial pattern and ecological analysis. Vegetatio 
80 (1989) 107-138. 

[22] Li X.-Y. Application of computation geometry in wireless networks. In: 
Cheng X., Huang X., Du D.-Z. (Eds.) Ad Hoc Wireless Networking (Kluwer 
Academic Publishers, 2004) 197-264. 

[23] Li X., Shi Y., Wang L., An updated survey on the Randic index. Math. Chem. 
Monogr. (6) (2008) 9-4 

[24] Magwene P. W. Using correlation proximity graphs to study phenotypic 
integration. Evolutionary Biology. 35 (2008) 191-198. 

[25] Matula D. W. and Sokal R. R. Properties of Gabriel graphs relevant to 
geographic variation research and clustering of points in the plane. Geogr. Anal. 
12 (1980) 205-222. 

[26] Muhammad R. B. A distributed graph algorithm for geometric routing in ad 
hoc wireless networks. J Networks 2 (2007) 49-57. 

[27] Plavsic D., Nikolic S., Trinajstic N. and Mihalic Z. On the Harary index for the 
characterization of chemical graphs. J. Math. Chem. 12 (1993) 235-250. 

[28] Randic, M. Characterization of molecular branching. Journal of the American 
Chemical Society 97 (1975) 6609-6615. 

[29] Runions A., Fuhrer M., Lane B., Federl P., Rolland-Lagan A.-G., and 
Prusinkiewicz P. Modeling and visualization of leaf venation patterns. ACM 
Transactions on Graphics 24 (2005) 702-711. 

[30] Sabljic A. and Protic M. Molecular connectivity: A novel method for 
prediction of bioconcent ration factor of hazardous chemicals. Chemico- 
Biological Interactions 42 (1982) 301-310. 



18 



[31] Santi P. Topology Control in Wireless Ad Hoc and Sensor Networks (Wiley, 
2005). 

[32] Sokal R. R. and Oden N. L. Spatial autocorrelation in biology 1. Methodology. 
Biological Journal of the Linnean Society 10 (2008) 199-228. 

[33] Song W.-Z., Wang Y., Li X.-Y. Localized algorithms for energy efficient 
topology in wireless ad hoc networks. In: Proc. MobiHoc 2004 (May 24-26, 
2004, Roppongi, Japan). 

[34] Sridharan M. and Ramasamy A. M. S. Gabriel graph of geomagnetic Sq 
variations. Acta Geophysica (2010) ,10.2478/sll600-010-0004-y, 

[35] Toroczkai Z. and Guclu H. Proximity networks and epidemics. Physica A 378 
(2007) 68. arXiv:physics/0701255vl, 

[36] Wan P.- J., Yi C.-W. On the longest edge of Gabriel Graphs in wireless ad hoc 
networks. IEEE Trans, on Parallel and Distributed Systems 18 (2007) 111-125. 

[37] Watanabe D. A study on analyzing the road network pattern using proximity 
graphs. J of the City Planning Institute of Japan 40 (2005) 133-138. 

[38] Watanabe D. Evaluating the configuration and the travel efficiency on proximity 
graphs as transportation networks. Forma 23 (2008) 81-87. 

[39] Yang Y., Lu L. The Randic index and the diameter of graphs Discrete 
Mathematics 311 (2011) 13331343. 



19 



